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The formation of interface dipoles in self-assembled monolayers (SAMs) of -CH3 and -CF3 ter- 
minated short-chain alkanethiolates on Ag(lll) is studied by means of density functional theory 
calculations. The interface dipoles are characterized by monitoring the change in the surface work 
function upon adsorption of the SAM. We compare results obtained for SAMs in structures with 
a different packing density of molecules, i.e. {^/l x \/7)i?19.1°, (\/3 x \/3)J?30°, and p(2x2). The 
work function of alkanethiolate SAMs on silver depends weakly on the packing density; that of 
fluorinated alkanethiolates shows a stronger dependance. The results are analyzed in terms of two 
nearly independent contributions to the interface dipole. These originate respectively from the 
molecular dipoles and from a charge transfer between the metal surface and the molecules. The 
charge transfer is determined by the silver-sulfur bond and it is independent of the electronegativity 
of the molecules. 



1. INTRODUCTION 

Self-assembled monolayers (SAMs) of organothiolates 
on metal surfaces are studied for a wide range of techno- 
logical applications running from catalysis, biosensors to 
microelectronic devices [T] [SJ [5] . In organic light-emitting 
diodes, the interfaces between the metal contacts and the 
organic material are critical in the device performance, 
since they control the injection of electrons and holes into 
the device |4]. Chemisorption of a SAM on a metal sur- 
face can alter its work function substantially. Depending 
on the SAM, the work function can be manipulated ad- 
vantageously to lower the energy barrier for electron and 
hole injection [3 |SJ |7] . 

Self-assembled monolayers have also become attractive 
for fundamental studies in metal-organic interfaces and 
molecular electronics. They represent stable and ordered 
structures, which can be prepared experimentally in air, 
in solution, or in vacuum [SJIS]. SAMs of alkanethiolates, 
C„H2n+iS, on Au(lll) are among the most extensively 
studied systems. Alkanethiolate SAMs on Au(lll) adopt 
a (-s/S X •\/3)-R30° structure or superstructures thereof 
[ini E] . Alkanethiolates form SAMs on a wide range of 
(noble) metal surfaces, which have a similar structure as 
on Au(lll). Variations on the packing density are possi- 
ble, however, and on Ag(lll) a somewhat denser packed 
(\/7 X \/7)-R19.1° structure has been reported [3]. 

The change in work function of the surface upon ad- 
sorption of the SAM is directly proportional to the dipole 
moment density generated at the SAM-metal interface. 
For SAMs on Au(lll) it has been shown that this dipole 
moment density is mainly determined by the permanent 
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dipoles in the thiolate molecular layer. The sulfur-gold 
bonds that are formed upon adsorption, are nearly apolar 
and give a very small contribution to the interface dipole 
[m [131 [HI [T5]. However, a small sulfur-metal bond 
dipole is typical of gold and the existence of a much larger 
bond dipole is indicated by experiments of alkanethiolate 
SAMs on silver [7] . In a previous computational study on 
model structures we have shown that large bond dipoles 
can be formed in the adsorption of SAMs on Ag and Pt 
surfaces [16]. 

In this paper we study the interface dipole formation 
resulting from adsorption of SAMs on Ag(lll) by first- 
principles density functional theory (DFT) calculations. 
In particular we examine the influence of the structure 
and the packing density of the molecules in the SAM. 
The (-\/7 X A/7)i?19.1° packing, which is observed exper- 
imentally for alkanethiolate SAMs on Ag(lll) [SI, is our 
starting point. We consider several low energy structures 
[n|[18]. The results are compared to the (Vs x \/3)i?30° 
structure, where the surface area per adsorbed molecule 
is 29% larger, which is the most common structure of 
alkanethiolate SAMs on other noble metal (111) surfaces. 
We also consider the less densely packed p(2x2) struc- 
ture, which has a 71% higher surface area per molecule. 
Fluorinated alkanethiolate SAMs on Au(lll) can be ob- 
served in this structure ^ , and it might be possible that 
this structure is also formed by such molecules adsorbed 
on Ag(lll). We show that although the interface dipole 
density is smaller for less densily packed structures, it 
is not simply proportional to the packing density due to 
dielectric screening in the molecular layer. 

The commonly used DFT functionals describe the for- 
mation of chemical bonds and the resulting charge dis- 
tribution very well, but they fail to capture the van der 
Waals interactions between the alkyl chains. Van der 
Waals interactions are relatively unimportant in short 
chain alkanethiolates, which is why we focus on the short 
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chain thiolates CH3S and C2H5S. To elucidate the influ- 
ence of the polarity of the molecules on the SAM— metal 
interface dipole, we also study the fluorinated thiolates 
CF3S and CF3CH2S. Since the directions of the dipole 
moment of fluorinated and of nonfluorinatcd thiolates are 
roughly opposite, this leads to an obvious difference in 
the interface dipole between SAMs of the two types of 
molecules. In addition, fluorinated thiolates have a much 
higher electronegativity. One would expect that by vary- 
ing the relative electronegativity of the surface and the 
molecules, one can modify the electron transfer between 
surface and molecules, which would give an additional 
contribution to the interface dipole. 

In this paper we will show that increasing the elec- 
tronegativity by fluorinating the alkyl tails does however 
not lead to a change in charge transfer. We will arrive at 
this conclusion by analyzing the interface dipole and sep- 
arating it into a contribution from the molecular dipoles 
and from the charge reordering at the metal— SAM inter- 
face. By comparing these results to those obtained for 
SAM-Au(lll) and SAM-Pt(lll) interfaces, it can be 
concluded that the charge transfer depends on the metal 
surface and the nature of the sulfur— metal bond, but not 
on the molecular tails. 

The paper is organized as follows. In the next sec- 
tion we describe the techniques we use for calculating 
and analyzing the interface dipoles and give details on 
the parameters used in the calculations. Subsequently 
the results on the SAM— Ag(lll) interfaces are discussed. 
First we discuss the possible structures and then we an- 
alyze the interface dipoles. The last section contains a 
short summary and the conclusions. 

2. THEORETICAL SECTION 

2.1. Total energy calculations 

The Ag(lll) metal surface is represented by a slab of 
layers of metal atoms stacked according to an fee ABC se- 
quence. A typical slab consists of four layers. The SAM 
is adsorbed on one side of the slab. The surface unit 
cell depends upon the monolayer structure and coverage. 
The cells used in our calculations are x ^/7)R19.1° , 
(VS X V3)i?30°, and p(2x2), which contain 7, 3 and 4 
metal atoms per layer, respectively. Periodic boundary 
conditions are applied in all three directions. This means 
that not only are the cells repeated along the surface 
(the xy-plane), but also the slabs are repeated in the z- 
direction. The atoms in neighboring cells are separated 
along the z-direction by a vacuum region of ~ 12 A. 
To cancel the artiflcial interaction between the dipoles of 
the repeated slabs, the Neugebauer-Scheffler dipole cor- 
rection is applied [19j . 

The electronic structure is treated within density func- 
tional theory (DFT) [20] using the PW91 functional [21] 
to describe the electronic exchange and correlation. The 
calculations are performed with the program VASP (Vi- 
enna ab initio simulation package) |22l [^j using the pro- 



jector augmented wave (PAW) method f^HI^S]. For no- 
ble metal atoms the outer shell s and d electrons are 
treated as valence electrons, and for first and second row 
elements the outer shell s and p electrons. The valence 
pseudo wave functions are expanded in a basis set con- 
sisting of plane waves. All plane waves up to a kinetic 
energy cutoff of 500 eV are included. 

The geometries are optimized by allowing the atoms in 
the top two metal layers and the atoms in the SAMs to 
relax. The (1x1) surface unit cell parameter is fixed at its 
optimized bulk value of 2.93 A. The calculations use a k- 
point sampling mesh of 7 x 7 for the {\/7 x y/7) structure, 
11 X 11 for the (VS x VS) and 9 x 9 for the (2 x 2) struc- 
tures, according to the Monkhorst-Pack scheme. For 
geometry optimization a Methfessel-Paxton smearing is 
used with a broadening parameter of 0.2 eV |2^. The 
energies of the optimized geometries are recalculated us- 
ing the tetrahedron scheme P7|. Tests regarding the 
slab thickness, vacuum thickness, k-point sampling grid, 
and plane wave kinetic energy cutoff are performed, from 
which we estimate that total energy differences are con- 
verged to within ~ 0.05 eV. 

The adsorption energy of the SAM is calculated by 
comparing the total energy of the slab (with the adsorbed 
SAM) , with that of the clean slab (with the top surface 
in its relaxed Ag(lll) structure), and the free alkanethi- 
olate (radical) molecules. If SAM adsorption results in a 
reconstruction of the surface that involves metal adatoms 
[TBI l28j , we assume that these adatoms are supplied by 
the bulk metal. The adsorption energy per molecule 
Ends associated with a surface structure that contains M 
molecules, TVs metal atoms per layer and A^ad adatoms is 
then given by 

-E'ads = -J^ [-E-slab ~ ^sEdcan — A^ad-E-bulk] — ^-mob (1) 

where E'siab is the total energy of the slab, -Ecican is the 
total energy of the clean slab per surface atom (top sur- 
face in optimized Ag(lll) structure), -Ebuik is the total 
energy of bulk Ag per atom, and E'moi is the total energy 
of an alkanethiolate molecule. Note that E'ads is negative 
if the adsorption is exothermic. 

To analyze the results we have also calculated sev- 
eral properties of isolated thiolate molecules: dipole 
moments, ionization potentials, electron affinities and 
electronegativities. For these calculations we use the 
GAMESS program and treat the electronic struc- 
ture within DFT using the B3LYP functional [501 131] ■ 
We use the aug-cc-pVTZ basis set. Calculations with 
the smaller 6-31 IG** basis set give dipole moments that 
are up to ~ 0.15 D smaller, and ionization potentials, 
electron affinities that differ by ~ 0.1 eV. 

2.2. Work functions and interface dipoles 

Interface dipoles can be extracted from the change in 
the surface work function upon adsorption of a SAM on 
a metal surface, as will be described below. Surface work 
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Figure 1. Plane averaged electrostatic potential V{z) 
of a slab consisting of four layers of silver atoms and 
one layer of CF3S molecules on top. The z axis is 
normal to the (111) surface. 



functions are evaluated from the expression: 

W ^Vioo) - Ef, 



(2) 



where V{oo) is the electrostatic potential in vacuum and 
Ep is the Fermi energy of the bulk metal, ^(oo) is ex- 
tracted by calculating the average electrostatic potential 
in the xy-planes of the slab: 



V{x, y, z)dxdy, 



(3) 



where C is the area of the surface unit cell and V{x, y, z) 
is the total electrostatic potential. The latter is gener- 
ated on an equidistant real space grid and the integral 
is obtained by straighforward numerical integration. In 
practice V{z) reaches an asymtotic value V{oo) within a 
distance of ^ 5 A from the surface into the vacuum. 

An example of V{z) is shown in Fig. [Ijfor CF3S on the 
Ag(lll) surface, l^metai is the work function of the clean 
Ag surface and VFsAM-Ag is the work function of surface 
covered by the SAM. Slab calculations produce a reason- 
able value for the bulk Fermi energy Ep, but a more accu- 
rate value is obtained from a separate bulk calculation. 
We follow the procedure described by Fall et al. [35]. 
From the convergence tests discussed in the previous sec- 
tion we estimate that calculated work functions are con- 
verged to within ~ 0.05 eV. Typically, DFT/PW91 cal- 
culations give work functions that are within 0.1-0.2 eV 
of the experimental values, although occasionally larger 
deviations of ^ 0.3 eV can be found [33| l34]. 

Upon adsorption of a SAM the work function of a metal 
surface usually changes considerably. The work function 
change AW can be interpreted in terms of a change in 
the surface dipole A/x: 



AW = WsAM- 



-metal 



metal 



(4) 



where A is the surface area per adsorbed molecule |35j . 
Note that Afj, corresponds to the component of the dipole 



moment directed along the surface normal, since only this 
component affects the work function. A/x is the result of 
the interface formation between the SAM and the metal 
surface, and we call it the interface dipole in the follow- 
ing. 

We split the interface dipole into a contribution usam 
from the molecular dipoles in the SAM, and a contribu- 
tion /ichom from the charge transfer between the metal 
surface and the molecule, which occurs upon chemisorp- 
tion of the SAM. The latter contribution is then defined 
by: 



(5) 



A'SAM is obtained from a separate calculation on a free- 
standing SAM without the presence of a metal slab, but 
with the molecules frozen in their adsorbed geometry. 
In the following we will show that /ichom is nearly in- 
dependent of the (fluorinated) alkyl tail of the thiolate 
molecule. This means that /ichom is mainly determined 
by the sulfur-metal bond and the charge transfer associ- 
ated with this bond. 

Note that the calculation of /isAM is done for a full 
monolayer. In practice /isAM is obtained from the ex- 
pression: 



/^SAM 



sqAAV 



(6) 



where AV^ = y(oo) — V{—oo) is the potential drop over 
the SAM, and V{oo),V{—(x:)) are the asymptotic elec- 
trostatic potentials on both sides of the SAM. These are 
easily obtained, since the potential reaches its asymptotic 
values within a distance of few A of the SAM. 

This calculation incorporates the effect of the depolar- 
izing electric field within the SAM that is generated by 
the close-packed molecular dipoles. Often this effect is 
modeled phenomenologically by introducing an effective 
dielectric constant e for the SAM: 



^J■SAM 



(7) 



where fiz is the z-component of the permanent dipole of 
the isolated molecule. By obtaining /i^ from a separate 
calculation we will extract the effective dielectric con- 
stant as a function of the packing density of the molecules 
in the SAM. 

. INFLUENCE OF PACKING DENSITY 

We will first discuss the possible structures of thiolate 
SAMs on the Ag(lll) surface and then study the influ- 
ence of the packing density on the work function. 

3.1. Structures 

From early scanning tunneling microscopy (STM) ex- 
periments it was concluded that SAMs of alkanethi- 
olates on the Ag(lll) surface form a commensurate 
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{V7 X V7)R10.9° structure [Ml |37]. Normal incident 
X-ray standing wave (NIXSW) experiments have con- 
firmed the \/7 X \/7 structure, but have corrected the 
registry of the SAM on the underlying substrate to 
{V7 X V7)Rl9.r [T7]. The proposed model of this 
structure has the molecules in the SAM arranged in 
a hexagonal lattice with a nearest neighbor distance 
between the sulfur atoms of 4.41 A, see Fig. |2ja). 
Long chain alkanethiolates adopt an expanded incom- 
mensurate {\/7 X v'7)^19.1° structure with a nearest 
neighbor distance of 4.6-4.8 A [38j [39], whereas short 
chain alkanethiolates keep the commensurate structure 
[Tfj . Recently a new model has been proposed for the 
(V7 X V7)R19.1° structure of CH3S SAMs on Ag(lll)on 
the basis of NIXSW and medium energy ion scatter- 
ing (MEIS) experiments [TH [101 [S] • It involves a sur- 
face reconstruction consisting of a 3/7 monolayer of Ag 
adatoms, which are bonded to the the methylthiolate 
molecules. Fig. [2j^c). 

The (a/? X \/7)Rl9.1° structure proposed first consists 
of three molecules per surface unit cell with two of the 
molecules adsorbed on hollow sites and one on a top site 
|17j . Starting from this structure we have relaxed the 
geometry of CH3S on Ag(lll) and the result is shown in 
Fig. [2]^ a) and (b). The molecules labeled 1 and 3 change 
their position only slightly and remain adsorbed in a hol- 
low site. The molecule labeled 2 moves away from the top 
site towards a bridge site. The angle between the surface 
normal and the C-S bond is 42°, whereas that angle for 
molecules 1 and 3 is only 9-10°. The latter molecules are 
almost standing upright, as can be observed in Fig. [2jb). 
We call this structure the "1,3 hollow" structure in the 
following. The most important bond distances and an- 
gles of this structure are given in Table [T] The adsorption 
energy (averaged, per molecule) according to Eq. ([l]) is 
-1.97 eV. 

In our previous calculations for alkanethiolate SAMs 
on Au(lll) the molecules show a strong preference for 
adsorption on bridge sites, instead of on hollow or top 
sites [TS]. Starting with CH3S molecules 1 and 3 on 
bridge positions we obtain the optimized geometry that is 
shown in Figs. [2je) and (f). In this structure the CH3S 
molecule 2 also moves closer to the bridge position, as 
compared to the 1,3 hollow structure, see Table [T] The 
new structure, which we call the "bridge" structure, is 
0.10 eV/molecule lower in energy than the 1,3 hollow 
structure. The bridge structure is energetically favored 
over the 1,3 hollow structure. The calculated adsorption 
energy in the bridge structure is —2.07 eV/molecule. The 
geometries of the three molecules in the bridge structure 
are more similar than in the 1,3 hollow structure. For 
instance, the angle between the surface normal and the 
C-S bond is in the range 43-51° for all three molecules, 
see Table [l| The work functions of the 1,3 hollow and 
the bridge structures are substantially different as will 
be discussed below. 

We have also optimized the geometry in the 
(\/7x V7)R19.1° reconstructed structure [TSlllI]. The 




Figure 2. Possible {V? x \/7)i?19.1° structures of 
CH3S SAMs on Ag(lll); (a), (b) top and side view of 
the 1,3 hollow structure; (c), (d) of the reconstructed 
structure; (e), (f) of the bridge structure; (g) of a 
C2H5S SAM in the bridge structure. "1, 2, 3" label the 
molecules, "a" labels the Ag adatoms. 



reconstruction involves a commensurate layer of Ag 
adatoms at a 3/7 monolayer coverage. The sulfur atoms 
of the adsorbed thiolate molecules are threefold coordi- 
nated by adatoms. The optimized geometry is shown in 
Figs, ^c) and (d) and bond distances and angles of this 
structure are given in Table [ij The resulting structure 
has the CH3S molecules standing upright with the C- 
S bond pointing along the surface normal in agreement 
with previous calculations [32] . The distance of the 
S atoms of the different molecules to the surface is the 
same within ~ 0.1 A, which represents a very small "rum- 
pling" in agreement with the latest experimental results 
[4T| . The work function of the reconstructed structure is 
substantially different from that of the unreconstructed 
structures as will be discussed below. The calculated ad- 
sorption energy in the reconstructed structure is —2.07 
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1,3 hollow 






molecule 1 


2 


3 


C-H(A) 


1.10 


1.10 


1.10 


S-C(A) 


1.83 


1.84 


1.84 


Ag-S(A) 


2.52/2.50/2.50 


2.47/3.01 


2.63/2.54/2.54 


C-S-normal(°) 


9.0 


42.0 


9.7 


Ag-S-Ag(°) 


78.7/77.5/77.4 


- 


76.2/71.3/72.5 


C-H(A) 




bridge 




1.10 


1.10 


1.10 


s-c(A) 


1.83 


1.84 


1.83 


Ag-S(A) 


2.56/2.49 


2.46/2.82 


2.58/2.48 


C-S-normal(°) 


51.4 


43.0 


46.6 


Ag-S-Ag(°) 


72.3 


- 


72.5 


C-H(A) 




reconstructed 




1.10 


1.10 


1.10 


s-c(A) 


1.84 


1.84 


1.84 


Ag-S(A) 


2.65 


2.66 


2.64 


C-S-normal(°) 











Ag-S-Ag(°) 


114.9 


115.3 


115.9 




2.47 


2.56 


2.47 



molecule 1 



TABLE I: Bond lengths and bond angles of the 
(V? X -\/7)i?19.1° 1,3 hollow, bridge and reconstructed 

structures of CH3S SAMs on Ag(lll). The columns 
indicate the three molecules in the supercell, see Fig. |2] 
dz is the distance along the surface normal between a 
Ag adatom and the top Ag layer. 



eV/molecule. This number is very close to the adsorption 
energy in the (unreconstructed) bridge structure. Within 
the intrinsic error bar of DFT calculations the two struc- 
tures are degenerate in energy. It has recently been sug- 
gested that the two structures, i.e. reconstructed and 
unreconstructed, might coexist on the surface [12]. 

The optimized unreconstructed {V7 x \/7)i?19.1° 
structure of a CH3CH2S SAM also has the molecules ad- 
sorbed on or near bridge sites, as shown in Fig. |2jg). 
The three molecules in the unit cell have a similar ge- 
ometry. For instance, the C-S-normal angle is 42-45°, 
see Table |TT] which is similar to the experimental value 
reported for ^3(^2)78 on Ag(lll) [IT]- The chain 
angles are in the range 12-17°, which is similar to exper- 
imental results for long-chain alkanethiolates . 

We have also studied (v^ x V3)i?30° and p(2x2) 
structures, where the surface area per adsorbed molecule 
is 29% and 71% larger, respectively, see Fig. [3| Experi- 
mentally it is not likely that alkanethiolates on Ag(lll) 
form these structures, but they enable us to model the 
influence of the packing density of the molecules in the 
SAM on the work function and the interface dipoles. In 
both these structure the bridge site is the favored ad- 
sorption site and the local geometry of the molecules is 
similar to that in the {^/7 x \/7)R19.1° bridge structure. 

The structure of (partially) fluorinated alkanethiolates 
is much less established than that of their nonfluorinated 
counterparts. On Au(lll) thiolates with long fluorinated 
alkyl tails have a less dense packing because of their rela- 
tively bulky tails For such SAMs a p(2x2) struc- 



C-H(A) 

c-c(A) 
s-c(A) 

Ag-S(A) 

C-C-normal(° 

C-S-normal(° 

chain(°) 

Ag-S-Ag(°) 



1.10 1.10 1.10 

1.52 1.53 1.52 

1.85 1.85 1.84 

2.57/2.50 2.47/2.90 2.58/2.49 



25.5 
44.2 
14.7 
72.1 



26.3 
41.9 
12.0 



26.8 
45.0 
17.1 
72.5 



TABLE II: Bond lengths and bond angles of the 
{V7 X V7)R19.1° bridge structure of CH3CH2S SAMs 
on Ag(lll). The chain angle represents the angle made 
by the line connecting the top C and S atoms with the 
surface normal. 



ture has been proposed, where the spacing between the 
adsorbate molecules is 5.87 A [IS]. SAMs of long-chain 
alkanethiolates with fluorinated end groups on Au(lll) 
have a (^3 x V3)m0° structure [IS]. The spacing be- 
tween the adsorbate molecules is then 5.08 A. 

We did not find reports on the structural details of 
fluorinated alkanethiolate SAMs on silver in the lit- 
erature. We optimized the structure of CF3S and 
CF3CH2S SAMs on Ag(lll) in three different packing 
densities, i.e., (^7 x V7), {VS x %/3) and p(2x2). For 
the {V7 X V7)R19.r structure of a CF3S SAM, which 
represents the most dense packing, we used the 1,3 hollow 
and the bridge structures of CH3S as starting points and 
optimized the geometry. The bridge structure of CF3S 
is more stable than the 1,3 hollow structure, albeit by 
less than 0.02 eV/molecule. However, even in the bridge 
structure only molecules 1 and 3 are actually adsorbed 
on bridge sites, whereas molecule 2 is adsorbed on a hol- 
low site. This results in a distorted hexagonal packing of 
the CF3S molecules. 

For the (^3 x V3)i?30° and p(2x2) structures of flu- 
orinated alkanethiolate SAMs on Au(lll) we have found 
a preference for the molecules to adsorb on bridge sites 
[TS]. The calculated nearest neighbor distance between 
the metal atoms on Au(lll) and Ag(lll) is very similar, 
i.e., 2.94 and 2.93 A respectively. Moreover, since even in 
the (\/7 X \/7) structure the molecules show a tendency 
to adsorb on bridge sites, we only consider bridge sites 
for the x ^3) and p(2x2) structures of SAMs on Ag, 
see Fig. |3] 

Table |III| lists the molecular geometries of CF3S and 
CF3CH2S SAMs on Ag(lll) adsorbed in bridge struc- 
tures for the different packings. The geometries are in 
fact very similar, with S-Au bond lengths in the range 
2.5-2.6 A, and angles of the C-S bond with the surface 
normal around 45°. In the case of CF3CII2S the chain an- 
gle is 13—15°, which is similar to that found in CH3CII2S, 
see Table [U 

The calculated adsorption energies indicate that a less 
dense packing of the fluorinated alkanethiolate SAMs is 
favorable, see Table III The p(2x2) structure is most sta- 



CFgS 

\/7 X V7 V3 X V3 p(2 X 2) 



C-F(A) 


1.35/1.36/1.35 


1.36 


1.36 


s-c(A) 


1.85/1.83/1.87 


1.84 


1.84 


Ag-S(A) 


2.51/2.41/2.57 


2.56 


2.54 


C-S-nomial(°) 


6.5/35.5/5.0 


43.2 


47.7 


Ag-S-Ag(°) 


79/-/72 


73.4 


71.5 


Sad. (eV) 


-1.98 


-2.32 


-2.39 


C-F(A) 




CF3CH2S 




1.35 


1.36 


1.36 


c-c(A) 


1.52 


1.52 


1.52 


C-H(A) 


1.10 


1.09 


1.09 


S-C(A) 


1.84 


1.85 


1.85 


Ag-S(A) 


2.53/2.45/2.52 


2.56 


2.56 


C-C-normal(°) 


25.9/25.7/29.0 


25.8 


21.0 


C-S-normal(°) 


44.0/42.4/42.7 


42.8 


46.7 


chain (°) 


13.9/12.9/15.1 


12.3 


16.7 


Ag-S-Ag(°) 


72.5/-/73.5 


72.1 


72.0 


Sad. (eV) 


-2.01 


-2.27 


-2.24 



TABLE III: Bond lengths, bond angles and 
adsorption energies of CF3S and CF3CH2S SAMs 
adsorbed on Ag(lll) surface in {\/7 x \/7)R19.1°, 
(\/3 X V3)m0° and p(2x2) structures. The chain angle 
represents the angle made in the CF3CH2S SAM by the 
top C atom and the sulfur atom with the surface 
normal. 



ble for CF3S, whereas for CF3CH2S the (^/3 x V3)R30° 
is slightly more stable. Some of the energy differences are 
quite small, but the calculations correspond to the exper- 
imental trend observed in fluorinated alkanethiolates on 
Au(lll) 051115]. 

3.2. Work functions and interface dipoles 

The change in work function upon adsorption of a SAM 
is defined by Eq. Q. For the clean Ag(lll) surface we 
calculate a work function of 4.50 eV, using a 25 x 25 
k-point Brillouin zone sampling grid. This is in good 
agreement with the experimental results of 4.5 eV |47| 
and 4.4 eV [7i], and with a previously reported theoretical 
value of 4.42 eV, which was extracted from DFT-GGA 
calculations and a 15 x 15 k-point sampling [48]. The 
calculated work function changes upon adsorption of the 



(fluorinated) alkanethiolate SAMs are given in Table IV 
The results clearly show that alkanethiolates decrease 
the work function, whereas fluorinated alkanethiolates in- 
crease it significantly. Depending on the structure, a dif- 
ference of 2-2.5 eV is found between the work functions of 
alkanethiolates and partially fiuorinated alkanethiolates 
adsorbed on the silver surface. The work function is thus 
tunable over a large range by adsorption of a suitable 
SAM, as is observed experimentally [SJ[7|. The absolute 
change in the work function upon adsorption of fiuori- 
nated alkanethiolates is 3-4 times larger than the change 
upon adsorption of nonfiuorinated molecules. This result 
is quite different from our previous findings for SAMs on 




Figure 3. CF3S SAM on Ag(lll) with the molecules 
adsorbed on bridge sites; (a),(b) top and side view of 
the (\/3 X V3)m0° structure; (c),(d) of the p(2x2) 
structure. 



Au(lll), where fluorinated and nonfluorinatcd molecules 
give a change in the work function that is similar in size 
(but of opposite sign, of course) p^. We conclude that 
the interaction between the molecules and the silver sur- 
face differs from that between the molecules and the gold 
surface. 

Kelvin probe measurements of the work function 
changes induced by adsorption of long-chain thiolate 
SAMs on Ag(lll) have been reported by Campbell et 
al. 5 , and by de Boer et al. [7], The values reported 
for CH3(CH2)9S and CH3(CH2)i5S are -0.7 eV [5], and 
-0.6 eV [7J, respectively, and 0.9 eV [S] and 1.1 eV [7] 
for CF3(CF2)7(CH2)2S. Our calculated value of -0.59 eV 
for CH3CH2S in the {V? x V7)R19.r bridge structure is 
close to the experimental values for alkanethiolates. The 
calculated values for the fluorinated molecules are higher 
than the experimental values even for the less densily 
packed p(2x2) structure. One explanantion might be 
that the effective dipole moment of a long-chain fluori- 
nated molecule embedded in a SAM is smaller than that 
of a short-chain molecule; in other words, the effective di- 
electric constant in a long-chain fluorinated thiolate SAM 
is larger, see Eq. ([T]). In addition, SAMs of fluorinated 
alkanethiolates may show more intrinsic disorder than 
their nonfluorinatcd counterparts [46 , which also reduces 
the average dipole moment perpendicular to the surface. 

The three x ^/7 structures for CH3S SAMs dis- 
cussed in the previous section give rise to different work 
functions. The adsorption of molecules in the bridge 
structure gives a substantial smaller shift of the work 
function than adsorption in the 1,3 hollow structure or 
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a/7xV7 V3xV3p(2x2) 

CH3S -0.52 (-0.84", -0.99") -0.61 -0.34 
CH3CH2S -0.59 -0.52 -0.46 

CF3S 1.80 (1.79") 1.61 1.48 

CF3CH2S 2.09 1.75 1.75 

TABLE IV: Work functions shifts in eV with respect 
to the clean surface of SAMs on Ag(lll) in the bridge 
structure; in the 1,3 hollow and the reconstructed 
structure, respectively. 

the reconstructed structure. Using Eq. ^ we can in- 
terprete the changes in the work function in terms of 
molecular dipoles. This difference in work function shift 
between the structures can be related to the orientation 
of the molecular dipoles. Whereas in the bridge struc- 
ture the molecular tails are tilted with respect to the 
surface normal, the tails of the two molecules that are 
adsorbed at hollow sites in the 1,3 hollow structure are 
almost perpendicular to the surface, see Table |l] The lat- 
ter leads to larger dipole moments along the surface nor- 
mal. All the molecular tails in the reconstructed struc- 
ture are perpendicular to the surface, which leads to large 
dipole moments and a large work function shift. These 
results suggest that work function measurements might 
be a simple experimental way of distinguishing between 
the different structures. 

As can be observed in Table |IV| the work functions 
have a relatively weak dependence on the packing den- 
sity. The local geometries of the thiolate molecules in the 
bridge structure are similar for the different packings, see 
e.g. Table Hence one would expect the individual 
molecular dipoles to be similar. The weak dependence 
of the work functions on the packing density might seem 
somewhat surprising. Assuming fixed interface dipoles 
Afi per molecule in Eq. (|4]), the work functions should 
scale as l/A, where A is the surface area per molecule. 
Since this is clearly not the case, it means that A/i de- 
pends on the packing density. In particular, the individ- 
ual molecular dipoles increase with decreasing packing 
density. 

Decreasing the packing density increases the distances 
between the molecules in the SAM. Hence it decreases 
the depolarizing field in the SAMs, or, in other words, 
the effective dielectric constant e introduced in Eq. ([7| 
decreases with decreasing packing density. This effec- 
tively increases the molecular dipoles, which opposes the 
effect of a decreasing density of the molecules on the in- 
terface dipole. The net result is a weak dependence of the 
work function on the packing density in the range con- 
sidered and in some cases even a nonmonotonic behavior, 
see Table HVl 

In order to quantify this analysis we make use of the 
relations given by Eqs. (|4])-([7]). We extract from the work 
function change an interface dipole per molecule A/i and 
split A/i into a contribution /is am from the molecular 
dipole and a contribution /ichcm from the charge transfer 



between the molecule and the surface upon chemisorp- 
tion. The results for each of the SAMs in the a/7 x -\/7, 
a/3 X a/3 and p(2x2) structures are reported in Ta- 
ble |Vj As we are explicitly interested in the influence 
of the packing density we continue to compare similar, 
i.e. bridge-like, structures. 

Ms AM is positive for fluorinated alkanethiolates, which 
means that the molecular dipoles point from the S atom 
to the CF3 group. For the nonfluorinated alkanethio- 
lates the absolute values of /isAM are larger, but the sign 
is negative, meaning that the dipoles point from the alkyl 
tails to the S atom. /<chom is positive for all molecules. 
The latter contribution is associated with a dipole that 
points from the surface to the molecule. It is associ- 
ated with a (partial) electron transfer from the surface 
to the molecule. Both contributions, /xsam and /ichcm, 
to the interface dipole Afj, are of comparable size. For 
the nonfluorinated molecules they are of opposite sign, 
which leads to moderate interface dipoles A/i ~ —0.2 to 
— 0.3 D and work function changes < 0.5 eV. The con- 
tributions /isAM and /ichcm have the same sign for the 
fluorinated molecules, which gives large interface dipoles 
A/i = 0.9-1.3D and large work function changes of up to 
2-2.5 eV. 

Comparing the results for the different packing den- 
sities in Table |V] we observe that the absolute value of 
fJ-SAM increases if the packing density decreases. This 
can be understood by noting that the effective dielec- 
tric constant e of the SAM decreases if the packing den- 
sity decreases, as discussed above. We have calculated 
the dipole moment of the isolated alkanethiolate radical 
molecules, fixing the molecules in the geometries they 
attain in the SAM. The component /i^ along the surface 
normal in the adsorbed geometry is given in Table |VI[ 
From Eq. ([t]) we then calculate the effective dielectric 
constant of the free standing SAM. The dielectric con- 
stants for fluorinated alkanethiolate SAMs are somewhat 
larger than those of their nonfluorinated counterparts. 
This might be expected since the polarizability of fluori- 
nated molecules is larger and therefore screening in the 
SAM is larger. The results in Table [V| clearly show that 
the dielectric constants decrease with decreasing packing 
density. 

Mchcm also increases with decreasing packing density. 
The origin of this effect is similar to that discussed in the 
previous paragraph; the screening of the dipoles in the 
layer decreases if the packing density decreases. An inter- 
esting observation is that at flxed packing density /ichcm 
shows little variation within the range of molecules. Ap- 
parently it is mainly determined by the S-Ag interaction 
and not so much by the molecular tails. 

Upon adsorption electronic charge is transferred from 
the surface to the molecule. In order to visualize the 
charge transfer at the interface upon adsorption of the 
SAMs, we calculate the change in electron density An: 

An = nsAM-Ag - ?^Ag - "-sam, (8) 
where nsAM-Ag, '^Ag and nsAM are the electron densities 
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V7 xV7 








X V3 






p(2x2) 






CH3S C2H5S CF3S CF3CH2S 


CH3S 


C2H5S 


CF3S 


CF3CH2S 


CH3S 


C2H5S CF3S CF3CH2S 


An 


-0.24 -0.27 0.83 


0.96 


-0.36 


-0.31 


0.95 


1.04 


-0.27 


-0.36 1.17 


1.38 




-0.78 -0.79 0.39 


0.43 


-0.88 


-0.79 


0.44 


0.50 


-0.92 


-1.01 0.50 


0.71 


/^chcm 


0.54 0.52 0.44 


0.53 


0.52 


0.46 


0.51 


0.54 


0.65 


0.65 0.67 


0.67 


e 


1.6 1.4 2.0 


2.7 


1.4 


1.4 


1.8 


2.3 


1.3 


1.1 1.6 


1.7 



TABLE V: Dipole per molecule Afj, from work function shift upon adsorption, the (perpendicular) molecular 
dipole moment /xsam in a free standing SAM and the chemisorption dipole moment ^chem of the SAMs on Ag(lll) 
surface. The values are in D. e is the effective dielectric constant of the free standing SAM. 



a) b) e) f) 




Figure 4. Difference electron density along the surface normal (z) averaged over the xy plane in unites of A 
and as isodensity surface, (a), (b) CH3S and (c), (d) and CF3S in the -\/7 x \/7 structure; (e), (f) CH3CH2S and (g), 

(h) CF3CH2S in the p(2x2) structure. 



of the SAM adsorbed on Ag(lll), of the Ag(lll) surface 
and of the free-standing SAM, respectively. The electron 
distributions are obtained on a real space grid from sepa- 
rate calculations on the adsorbed SAM, on the substrate 
and on the free-standing SAM, respectively. In the lat- 
ter two calculations the substrate and the molecules are 
fixed in the adsorption geometries. 

As examples. Fig. |4] shows the difference electron den- 
sity An averaged in the xy plane along the surface normal 
of CH3S and CF3S SAMs adsorbed in the ^7 x \/7 struc- 
ture and of CH3CH2S and CF3CH2S SAMs adsorbed in 
the p(2x2) structure. In addition three dimensional vi- 
sualizations of Art at the interface are presented. Only 
the region around the SAM/Ag(lll) interface is shown, 
since in the substrate and the vacuum region An ~ 0. 
The figures clearly demonstrate that An is localized at 
the interface, i.e., near the sulfur atoms and the top metal 
layer. Electronic density is transferred mainly from the 
top layer of silver atoms to the sulfur atoms, which results 
in a dipole moment /Uchem- 

To check the consistency of this analysis, we can calcu- 
late the metal-sulfur dipoles from the difference electron 



density: 



A^chcm 



2An(x, y, z)dxdydz, 



(9) 



where we choose Zg in the center between the second and 
the third metal layer of the substrate and Zy in the center 
of the vacuum. The metal-sulfur dipoles fJ-'rh^,,^ are within 
10% of the values /ichcm listed in Table [Vj 

By integrating the peak of An on the sulfur atom, see 
Fig. |4] one can calculate the charge transfer from the 
substrate to the sulfur atom. A typical value over a range 
of structures is g = (0.11 ± 0.01)e. Modeling the charge 
transfer dipole as Hchem — Q_d gives d — 1.1 A, assuming 
A'chcm — 0.6 D. The distance between the sulfur atoms 
and the top layer of silver atoms is 2.0-2.2 A, so this 
analysis is consistent with the interpretation of /ichcm as 
a metal-sulfur bond dipole. 

The fact that this dipole moment hardly depends on 
the different molecular tails is slightly surprising, since 
the electronegativity of fluorinated tails is much higher 
than that of unfluorinated ones. The MuUiken elec- 
tronegativity of a molecule is defined as: 



EA + IP 



XM 



(10) 



9 



CHsS C2H5S CF3S CF3CH2S 


EA 


1.73 


1.87 3.02 


2.43 


IP 


9.20 


8.95 10.79 


9.82 


Xm 


5.47 


5.41 6.91 


6.13 




1.70 


1.82 1.05 


2.06 




-1.23 


-1.14 0.79 


1.17 



TABLE VI: Electron affinity (EA), ionization 
potential (IP), MuUiken electronegativity (xm) in cV, 
total dipole moment |/itot| and dipole moment along the 
surface normal /i^ in D of isolated molecules in their 
adsorbed geometries. 

where EA and IP are the electron affinity and ionization 
potential of the molecule. The EAs, IPs and MuUiken 
electronegativities of the molecules considered in this pa- 
per (fixed in their adsorbed geometries) are given in Ta- 
ble |VT1 

One observes a considerable difference in the elec- 
tronegativities Xm of the molecules. The xai of both 
alkanethiolates is similar, but the xm of fluorinated alka- 
nethiolates is much larger. The HOMO of the radical 
neutral molecules, which plays a role in determining the 
EA and IP, is stabilized by the electron withdrawing 
CF3 group. This property is commonly associated with 
the attractive Coulomb field of the CF3 group. The 
HOMO of the neutral molecules is localized mainly on 
the sulfur atom and one expects that the effect of the 
CF3 group decreases if the distance between this group 
and the sulfur atom increases. Indeed one finds that the 
EA, IP and xm of CF3S are significantly higher that 
those of CF3CH2S. 

The electronegativity of a metal surface is given by 
its work function. From simple chemical reasoning 
one would assume that the charge transfer between a 
molecule and a surface would depend on the difference of 
their electronegativities. This is clearly not the case; the 
electronegativity of the molecule does not seem to infiu- 
ence the charge transfer. This suggests that the effects 
of the Coulomb field of the CF3 group and the alkyl tails 
on the charge distribution at the sulfur-metal interface 
are screened by the metal. 

4. SUMMARY AND CONCLUSIONS 

We have studied the interface dipole formation and 
work function changes produced by adsorption of CH3S, 
CH3CH2S, CF3S and CF3CH2S SAMs on the Ag(lll) 
surface by means of DFT calculations. Adsorption of the 
alkanethiolates CH3S and CH3CH2S decreases the work 
function as compared to the clean metal surface, whereas 
adsorption of the fiuorinated alkanethiolates CF3S and 
CF3CH2S increases the work function. 

In particular we have examined the infiuence of the 
structure and the packing density of the molecules in 
the SAM. CH3S on Ag(lll) in the unreconstructed 
(\/7 X y/7)Rl9A° structure with two of the three 



molecules in the unit cell adsorbed at a hollow site, leads 
to a work function shift of —0.8 eV. Adsorbing the CH3S 
molecules on bridge sites stabihzes the structure by 0.10 
eV/molecule and gives a work function shift of —0.5 eV. 
The recently proposed surface reconstruction induced by 
CH3S adsorption yields an almost identical adsorption 
energy, and a work function shift of —1.0 eV. The differ- 
ence between the work functions of these structures can 
be interpreted in terms of the difference in the orientation 
of the molecular dipoles. 

These results are compared to the less densily packed 
{y/3 X V3)i?30° and p(2x2) structures, which are more 
likely to occur for fluorinated alkanethiolate SAMs. Al- 
though the work function shift generally decreases for de- 
creasing packing density, it is not simply proportional to 
the density of molecular dipoles. A partial compensating 
effect is caused by a decrease of the dielectric screening 
in the molecular layer. 

Comparing the different molecules adsorbed in similar 
geometries shows that fluorinated alkanethiolates can in- 
crease the work function by up to 2 eV, which is much 
larger that the decrease in work function caused by (non- 
fluorinated) alkanethiolate adsorption. We explain this 
by separating the interface dipole into a contribution 
from the molecular dipoles and from the charge reorder- 
ing at the metal— SAM interface. Electron transfer oc- 
curs from the Ag surface to the sulfur atoms of the 
thiolate molecules. The resulting dipole points in the 
same direction as the molecular dipole for fluorinated 
molecules. Addition of the two dipoles leads to a large in- 
terface dipole and a large work function shift. The direc- 
tion of the molecular dipole of nonfluorinated molecules 
is opposite to the metal-sulfur bond dipole, resulting in 
a much smaller interface dipole and work function shift. 

The electron transfer from the Ag surface to the 
molecules is remarkably independent of the electroneg- 
ativity of the molecules. In good approximation the 
charge reordering only depends upon the metal-sulfur 
bond, which suggest that the influence of the molecu- 
lar tails is screened by the metal substrate. In previous 
calculations we arrived at the same conclusion for adsorp- 
tion of alkanethiolate SAMs on other noble metal sur- 
faces, indicating that this result is more general [T^ITB]. 
For adsorption on Ag(lll) we flnd effective Ag-S dipoles 
^^hom = 0.51 ±0.04 D and 0.66±0.01 D in the (^/3 x ^) 
and p(2x2) structures, respectively. 

For adsorption on Au(lll) we have found very small 
Au-S dipoles ^chcm < 0.1 D, indicating an apolar Au- 
S bond [15 , whereas for adsorption on Pt(lll) the Pt-S 
dipole is /ichom = — 0.45±0.03 D The latter indicates 
an electron transfer from the sulfur atoms to the Pt sur- 
face. The metal-sulfur bonds formed upon SAM adsorp- 
tion generate dipole moments that reduce the work func- 
tion differences between the clean metal surfaces, and can 
even reverse the order. The molecular dipole moments 
exhibit similar values for adsorption on Ag, Au and Pt, 
thus giving the possibility to design interface dipoles. By 
adding molecular and metal-sulfur bond dipoles the over- 
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all work function can be determined. It is possible to ma- 
nipulate metal work functions considerably using SAMs. 
Work function shifts that can be as large as 2 eV. 
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